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ABSTRACT 

We present the first systematic investigation of the morphological and timing properties of flares in 
GRBs observed by Swift/XRT. We consider a large sample drawn from all GRBs detected by Swift, 
INTEGRAL and HETE-2 prior to 2006 Jan 3 1 , which had an XRT follow-up and which showed significant 
flaring. Our sample of 33 GRBs includes long and short, at low and high redshift, and a total of 69 flares. 
The strongest flares occur in the early phases, with a clear anti-correlation between the flare peak intensity 
and the flare time of occurrence. Fitting each X-ray flare with a Gaussian model, we find that the mean 
ratio of the width and peak time is (At/t) = 0.13 + 0.10, albeit with a large scatter. Late flares at times 
> 2000 seconds have long durations, At > 300 s, and can be very energetic compared to the underlying 
continuum. We further investigated if there is a clear link between the number of pulses detected in the 
prompt phase by BAT and the number of X-ray flares detected by XRT, finding no correlation. However, 
we find that the distribution of intensity ratios between successive BAT prompt pulses and that between 
successive XRT flares is the same, an indication of a common origin for gamma-ray pulses and X-ray 
flares. All evidence indicates that flares are indeed related to the workings of the central engine and, 
within the standard fireball scenario, originate from internal shocks rather than external shocks. While all 
flares can be explained by long-lasting engine activity, 29/69 flares may also be explained by refreshed 
shocks. However, 10 can only be explained by prolonged activity of the central engine. 

Subject headings: gamma rays: bursts — X-rays: bursts 



Introduction 



The advent of Swift (Geh rels et all 2004) has 



brought substantial advances in our knowledge of 
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GRBs including the discovery of the first afterglow 
(with a position known to several arcsec precision) 
of a short burst. Swift also brought on the defini- 



tion o f a possible third class of GRBs (IGehrels et al. 
2006), the discovery of a smooth transition between 



prompt and afterglow emission ( TagliaferrLetalj2005| 



IVaughan et al.ll2006HO'Brien et aU2006l). and the def- 
inition of a canonical X-rav light curve (jNousek et al. 



20061: lO'Brien et al.ll2006tlzhang et alJl2006l) . The lat- 
ter includes a steep early part (oc r a ' with 3 < a\ < 5, 
typically interpreted as GRB high-latitude emission), 
a flat phase (0.5 < a% < 1, generally interpreted as 
due to energy injection into the external shock), and 
a last, steeper part (1 < 03 < 1.5, the only one ob- 
served by pre-Swift X-ray i nstruments), wit h the pre- 



dicted t decay (see, also IWu et alj 120051) - Some 



times, a further steepening is detected after the nor- 



1 



mal decay phase, w hich is consistent with a jet break 
dZhang et al.ll2006h . 



What may be the most surprising discovery is 
the presence of flares in a large percentage of X-ray 
light curves. Fl ares had been p reviously observed 
in GRB 97050 8 dPiro et all 1 199% . GRB 011121 and 



GRB 011211 dPiro et alJ I2005IT IPiroetal.1 d2005l) 



suggested that the X-ray flares observed in the lat- 
ter two events were due to the onset of the afterglow, 
since the spectral parameters of these flares were con- 



2006b), GRB 



(Burrows etaH l2005bl 


Romano et al. 


050502B (] 


-alconeetalj |2006j), and 


(Pagani et al. 


20061), we 


have learned 



that flares can be considerably energetic and that they 
are often characterized by large flux variations. In- 
deed, the flare fluences can be up to 100% of the 
prompt fluence and the flare fluxes, measured with re- 
spect to the underlying continuum, AF/F, can vary 
in very short timescales Af/f pea k (AF/F ~ 6, 500 
and 25, Af/f peak « 1, ~ 1, ~ 1 in XRF 050406, 
GRB 050502B and GRB 050607, respectively, where 
Af measures the duration of the flare and f pea k is 
measured with respect to the trigger time). Further- 
more, detailed spectral analysis has proven that these 
flares ar e spectrally harder th a n the underlying con- 
tinuum dBurrows et alj l2005bl ; iRomano et alj |2006b; 



Falcone et al. 2006). In particular, they follow a hard- 



to-soft evoluti on, which is remi niscent of the prompt 
emission (e.g. iFord et alJll995l). The spect ra of the 



flares in GRB 050502B dFalcone et alj|2006l) are bet 



ter fit by a Band function dBand et alj 1993 , which 



is the standard fitting model for GRB prompt emis- 
sion), than by an absorbed power law (which usu- 
ally suffices for a standard afterglow). Very often 
multiple flares are observed in the same light curve, 
with an underlying afterglow consistent with having 
the same slop e before and after the flare. Finally, 
GRB 050724 dBarfhelmv et al.ll2005A [Campana et all 
2006h and GRB 050904 dCusumano et alJl2006h have 
demonstrated that flares happen in short GRBs as well 
as long ones, at low and very high redshift (the record 
being held by GRB 050904 at z = 6.29). 

The picture that the early detections o f flares have 
drawn was described by iBurrows et al. (2006) and 
Chincarini et al. d2006i and references therein), and 
a few conclusions were derived, albeit based on a 
small sample of objects. The presence of an under- 
lying continuum consistent with the same slope be- 
fore and after the flare (GRB 050406, GRB 050502B) 



seems to rule out external shocks models, since no 
trace of an energy injection can be found; the large 
observed AF/F cannot be produced by synchrotron 
self-Compton in the reverse shocks; the very short 
timescales Af/? pea k < 1 also generally rule out exter- 
nal shocks, unl ess very carefully b alanced conditions 
are met (e.g., iKobayashi et all2007h : furthermore, the 
flare spectral properties (harder than the underlying af- 
terglow, evolving from hard to soft) indicate a different 
physical mechanism from the afterglow, and possibly 
the same as the prompt one. 

In this work we present the first comprehensive 
temporal analysi s of all GRBs observe d by the X-ray 
Telescope (XRT, Burrows et al. 2005a)-both long and 
short, independently of whether they are GRBs, X- 
ray R ich (XRR) or X-ray Flashes (XRF, iHeise et al. 



2001) at low and high redshift-that showed flares in 
their X-ray light curves. We assess whether the ev- 
idence for prolonged engine activity accumulated on 
the first observed flares survives statistical investiga- 
tion and discuss the case that flares are indeed related 
to the workings of the central engine. We also present 
the results of a cross-check analysis between X-ray 
flares and pu lses detected by the B urst Alert Tele 
scope (BAT, iBarfhelmv et al. 2005 eh i n the gamma 



ray p rompt emission. A second paper dFalcone et al. 



2007) will study the same sample from the spectro- 
scopic point of view, in a natural complement to this 
work. 

This paper is organized as follows. In $2] we de- 
scribe our GRB sample and in §[3] the data reduction 
procedure; in $4] we describe our XRT data analy- 
sis and in $5] our cross-check analysis between X- 
ray flares and pulses detected by BAT in the gamma- 
ray prompt emission. In $6] we present our main re- 
sults and in $7] we discuss our findings. Throughout 
this paper the quoted uncertainties are given at 90% 
confidence level for one interesting parameter (i.e., 
Ax 2 = 2.71) unless otherwise stated. Times are re- 
ferred to the BAT trigger Tq, t = T — To, unless oth- 
erwise specified. The decay and spectral indices are 
parameterized as F(v, f) oc t~ a v~P, where F v (erg cm -2 
s' 1 Hz -1 ) is the monochromatic flux as a function of 
time t and frequency v; we also use Y — f3 + 1 as the 
photon index, N(E) oc E~ r (ph keV~' cirT 2 s ). 

2. Sample definition 

We considered all GRBs detected by Swift, INTE- 
GRAL and HETE-2, between the Swift launch and 
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2006 January 31 (119 events) for which XRT obtained 
a position (99). We then examined all light curves 
searching for large scale activity in excess of the un- 
derlying power-law light curve (flares). We defer a 
detailed analysis of small scale and small frequency 
variability, sometimes referred to as "flickering" to a 
later paper. None of the INTEGRAL- or HETE-2- 
triggered bursts showed any flares although we note 
that these bursts were observed by XRT much later 
than the Swift-triggered ones. As will be discussed in 
34.41 where we investigate the sample biases in depth, 
we evaluate the completeness of our sample with a 
large set of simulations. We established that our flare 
sample can be considered complete with respect to 
faint flares only at late times (typically 10 3 seconds 
after the trigger). In Table Q] we list all the GRBs that 
were selected for the analysis along with their redshifts 
(when available, i.e., for 9 of them), T^s, and BAT fiu- 
ences. This is what we shall refer to as our "full" sam- 
ple, consisting of 33 GRBs, on which we attempted 
the timing analysis described in §|4] The light curves 
of the full sample are shown in Fig.Q] 

Some light curves, however, were not fit for the full 
analysis. For instance, although joi nt analysis of BAT 
and XRT data on GRB 0502 19A dGoadet al.l l2006b 
showed a simultaneous flare (hence its inclusion in our 
sample), the portion of the flare that was observed with 
XRT was not long enough to fully characterize it. In 
the same manner, a handful of events (GRB 050826, 
GRB 051016B, GRB 060109) which are included in 
our full sample because they showed either low-signal 
late-time flares or a flattening in the XRT light curve, 
were excluded from a full analysis because of the low 
statistics obtained. All these special cases are reported 
in Table[T]in italics. After these exclusions, we defined 
our "restricted" sample, which consists of 30 GRBs on 
which we succeeded in performing our full analysis. 

We not e that our restricted sample differs from 
the one of iFalcone et al.l (120071) . because of different 
requirements for the analys is. As an example, for 



GRB 050820A IFalcone et al.l (120071) could perform de- 
tailed spectroscopic analysis of the flare portion ob- 
served by XRT, but our full timing analysis was not 
applicable. 

3. Data Reduction 

The XRT data were first processed by the Swift 
Data Center at NASA/GSFC into Level 1 products 
(event lists). Then they were further processed with 



the XRTDAS (vl.7.1) software package, written by 
the ASI Science Data Center (ASDC) and distributed 
within FTOOLS to produce the final cleaned event 
lists. In particular, we ran the task xrtpipeline 
(vO.9.9) applying calibration and standard filtering and 
screening criteria. An on-board event threshold of 
-0.2 keV was applied to the central pixel of each event, 
which has been proven to reduce most of the back- 
ground due to either the bright Earth limb or the CCD 
dark current (which depends on the CCD temperature). 

The GRBs in our sample were observed with differ- 
ent modes, which were automatically chosen, depend- 
ing o n source count rates, to minimize pile-up in the 
data et al. 2004). For the GRBs observed during 
the calibration phase, however, the data were mainly 
collected in Photon Counting (PC) mode, and pile-up 
was often present in the early data. Furthermore, for a 
few, especially bright GRBs (which were observed af- 
ter the Photo-Diode (PD) mode was discontinued due 
to a micrometeorite hit on the CCD) the Windowed 
Time (WT) data were piled-up, as well. Generally, 
WT data were extracted in a rectangular 40x20-pixel 
region centered on the GRB (source region), unless 
pile-up was present. To account for this effect, the WT 
data were extracted in a rectangular 40x20-pixel re- 
gion with a region excluded from its centre. The size 
of the exclusion region was determined follow ing the 
procedures illustrated in lRomano et al. d2006al) . To ac- 
count for the background, WT events were also ex- 
tracted within a rectangular box (40x20 pixels) far 
from background sources. 

The PC data were generally extracted from a cir- 
cular region with a 30-pixel radius. Exceptions were 
made for bright sources, which required a > 30-pixel 
radius, and for faint sources, which required a smaller 
radius in order to maintain a high signal-to-noise ratio. 
When the PC data suffered from pile-up, we extracted 
the source events in an annulus with a 30-pixel outer 
radius and an inner radius, depending on the degree of 
pile-up as determined via the P SF fitting method illus- 
trated in IVaughan et al. I d2006l) . PC background data 
were also extracted in a source-free circular region. 

For our analysis we selected XRT grades 0-12 and 
0-2 for PC and WT data, respectively (ac cording to 



Swift nomenclature; iBurrows et al. 2005a). To cal 



culate the PSF losses, ancillary response files were 
generated with the task xrtmkarf within FTOOLS, 
and account for different extraction regions and PSF 
corrections. We used the latest spectral redistribution 
matrices in the Calibration Database maintained by 
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HEASARC. 

From both WT and PC data, light curves were cre- 
ated in the 0.2-10keV energy band using a criterion of 
a minimum of 20 source counts per bin, and a dynam- 
ical subtraction of the background. Therefore, in our 
sample, each light curve was background-subtracted, 
and corrected for pile-up, vignetting, exposure, and 
PSF losses. 

4. Data Analysis 

The first goal of this work was to obtain a quan- 
titative assessment of flare characteristics. We thus 
set to measure statistical parameters such as the ratio 
of the flare duration to the time of occurrence Af/f, 
the power-law decay slope a^w, the decay to rise ratio 
A?f a ii/A?,i se , the flare energetics, and the flare to burst 
flux ratio. Different approaches suited the data best, 
depending on the flare statistics, as we outline below. 

4.1. Equivalent widths 

We calculated the equivalent width of the flares 
defined as EW = f dt, where 

« continuum (/ ) 

^continuum(f) describes the assumed shape of the contin- 
uum light curve underneath the flare (the local power 
law "underlying continuum") and F O b se rved(0 is the 
observed light curve, i.e., the combination of the con- 
tinuum and flare (the analytical fits to the continua 
are described in detail in 34. 21 and their parameters re- 
ported in Table 12). The equivalent width (expressed 
in units of seconds, as reported in Table [3] column 6) 
represents the time needed for integration of the con- 
tinuum to collect the same fluence as of the flare and 
it can give us a first indication of the lowest fluence 
we are able to measure for a flare. Indeed, the faintest 
equivalent width measured, on a rather weak afterglow 
with XRT fluence of ~ 1 .3 x 10~ 8 erg cirT 2 light curve, 
is 7.9 s in a small flare detected in GRB 050819. At 
the other extreme of the EWs is GRB 050502B, where 
we detect two flares, both characterized by large EWs. 
The first one is extremely bright and indeed has a 
fluence that is larger than the fluence of the underly- 
ing continuum light curve (1.43 x 10~ 6 erg cirT 2 and 
1.23 x 10~ 6 erg cirT 2 , respectively). Even though (see 
{34.4b our completeness for faint flares is somewhat 
limited at early times, this may be an indication that 
the flare is generally stronger than the continuum light 
curve and possibly an unrelated phenomenon. 

Our ability to measure EWs is limited by the dis- 
crete sampling of the light curves as well as the relative 



faintness of the flares, therefore we could only obtain 
EWs for 48 flares. Figure [2] shows the distribution of 
the EWs for our sample. 

4.2. Af/f from Gaussian fits 

The simplest analytical characterization of the flare 
morphology is obtained by adopting a multiply-broken 
power law to model the underlying continuum, and 
a number of Gaussians to model the superimposed 
flares. We adopted the following laws for the con- 
tinuum, i) simple power law: F(t) - Kr a> , ii) bro- 
ken power law: F(t) = Kt~°" for f < t bl and F(t) = 
Kt b "' (t/t bl y a2 for f > fbi, in) doubly-broken: F(t) = 
Kr ai for f < t M and F(t) = Kt~" 1 (f/fbiT" 2 for 
t M <t < f b2 , F(f) = Kt b ^(t b2 /t bl r a Ut/t b2 r a ' for 
f > fb2, and so on, where fbi and fb2 are the times 
of the breaks. For our flares, we iteratively added as 
many Gaussians as required to accommodate the x 2 
locally around each flare. The best-fit model parame- 
ters for each component (continuum and flares) were 
derived with a joint fit and are reported in Table [2] 
(continuum parameters) and Table [3] (flare parameters, 
columns 2-4). Column 5 of Table [3] reports flare peak 
fluxes measured with respect to the underlying contin- 
uum, or AF/F). The full gallery of fits is illustrated in 
Fig. Q] In Fig. [3] we show the distribution of the peak 
times (i.e., the Gaussian peaks). 

Based on these fits, we calculated Af/f for each 
flare, adopting the Gaussian width (cr) and peak t pes k 
as Af and f, respectively, where f pea k ranges between 
95 s and ~ 75 ks. We do not include the Gaussian 
fits for GRB 060124 for sake of homogeneity of the 
sample, since the XRT data include the prompt phase 
dRomano et alJl2006ah . Our ability to fit flares with 
Gaussians is less affected by discrete sampling of the 
light curves than for the EW determination, but it still 
suffers from the faintness of the flares; therefore we 
obtained fits for 69 Gaussian-modeled flares. In Fig. [4] 
we show the distribution of the Af/f, which peaks 
at 0.13, and which yields a mean value of (Af/f) = 
0.13 + 0.10. An assessment of selection effect that may 
affect this result is reported in 94.41 

4.3. Decay slopes, rising and decaying times from 
more realistic models 

Flare profiles can be quite complex. As an example, 
in Fig. [5] we show the light curves of GRB 050730, 
in which different flares are best fit by different laws 
(two power laws for the first, and an exponential rise 
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followed by a power-law decay for the second one), 
of GRB 050502B (first flare), and GRB 0601 11A. A 
more realistic description of the flare profile should 
therefore account for the skewness observed in many 
flares as well as different rising and falling slopes and 
times, which we indicate with a r [ se , omi, A?ri se , and 
A?f a u, respectively. Such a fit can be performed with 
power laws {F(t) oc (t - to)~ a ), in which case it is crit- 
ical to define the reference time to. In practice, for ?o 
we consider the peak time as well as the times, before 
and after the peak, when the flare profile deviates sig- 
nificantly from the continuum fit: Af^e = f pe ak - h, 
Affaii = ?2 - fpeak, where fi and ?2 are the times when 
the flare emission (fitted in this case by a Gaussian) 
crosses the fraction / of the flare peak emission on ei- 
ther side of the peak. For the calculation of the decay 
slopes, we chose / = 0.01. As for the previous fits 
( 14.11 14.21 ). our procedure requires a power law contin- 
uum beneath each flare. The values of am we com- 
puted for this sample (consisting of 35 flares) are re- 
ported in Table [3] (column 7) and their distribution is 
shown in Fig. [6] We derive (artyi) = 3.54 with stan- 
dard deviation of <x = 1 .50. We note that our choice 
of / = 0.01 was an operative decision; using a differ- 
ent definition for t\ and t%, the measure of the slope 
decay changes. For instance, for the large flare in 
GRB 050502B we obtain a m = 6.32, 5.58, and 5.20 
for / = 0.01, 0.05 and 0.10, respectively. This shows 
how critical the definition of t\ and ?2 is in measuring 
the decay slope that describes the temporal behaviour 
of a burst or flare. 

The quantities Af r ; se and Affaii are in themselves 
qu ite interesting , since , as is well known from the work 
of iNorris et al.N 19961) a nd fr om the simulations by 
Daigne & Mochkovitchl (119981) . the observed bursts, 
which are due to the internal shocks, present a Fast- 
Rise Exponential-Decay (FRED) shape with a ratio 
Aff a ii / Afri se = 3.4. We calculated Afrj se and Aff a u for 
our sample defining them in terms of / = 0.05, by as- 
suming the underlying power law continuum beneath 
the flares, and performing a separate fit to the rising 
and decaying part of the flare light curve. While for 
the decaying part we always used a power law, we 
found that in many instances the best fit to the rising 
part was obtained with an exponential. Using these 
fits, we calculated T90 (the time defined by / = 0.05) 
and the ratio At ft adopting Af = T90 and t - f pea k- 
Table |3] reports T90, Affaii /Afrise, At/t (columns 8-10), 
while Fig.|7]shows the distributions of Aff a ii/Af r ; se . 



4.4. Selection effects 

As stated above, our ability to measure statistical 
quantities from the light curves critically depends on 
both the discrete sampling of the light curves as well 
as the actual intensity of the flares with respect to the 
continuum beneath them. In this section we present 
our considerations on the biases that may affect our 
analysis and their effect on our results. One of the first 
difficulties comes from the blending of flares, which 
causes the EW, At, At/t, to be overestimated. Our re- 
sult of low At/t is thus an upper limit on the intrinsic 
sharpness of flares. 

4.4.1. Time resolution and low-earth orbit biases 

The time resolution of our observations, which 
decreases logarithmically during the XRT afterglow 
follow-up, is the first critical factor. Typically, at the 
beginning of the XRT light curve the sampling is quite 
good, but if the flare duration is of the order of the time 
it takes it to fade, then it will not be possible to rec- 
ognize it as such, and it will be interpreted as a steep 
power law, instead. This was often observed in the 



early X RT light curves, as reported by Tagliaferr i et al. 



(12005b and lO'Brien et alj (120061) . and it is partially re- 
lated to the short but finite time (usually > 60 s) it takes 
Swift to re-point to the GRB. On the other hand, at the 
end of the XRT light curve, the sampling also degrades 
because of the long integration required to achieve suf- 
ficient S IN, so that flares shorter than the integration 
time are smeared out and consequently, except for the 
brightest ones, their resulting average count rate drops 
below the detection threshold. 

Due to Swift's low-earth orbit, the data are not col- 
lected in a continuous way, but in portions of an orbit 
that last less than an hour. This is illustrated in Fig. [8] 
(left), which represents the distribution of the observ- 
ing times relative to the BAT trigger, of all the light 
curves in our sample. For each observation of the light 
curve, we estimated the time, which we shall refer to 
as bin time (BT), within which the counts were accu- 
mulated in order to have a S IN > 3. For t > 10 4 s 
the BT will generally include data from consecutive 
orbits. In Fig. [8] (right), we show the time resolution 
(BT) as a function of the time since the BAT trigger, as 
well as the curve that corresponds to BT/f = 0.1 and 
lies above the large majority of the data. It indicates 
that the instrumental resolution BT/t = (At/t)s/N=3 is 
in most cases significantly better than At/t sa 0.1 and 
is often even better than 0.01. In other words, our data 
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are not biased against At/t < 0.1. 

4.4.2. Biases in the sample definition criteria 

In order to evaluate the completeness of our sample 
we tested the sample definition criteria against selec- 
tion effects by means of simulations. 

First of all, for each flare in our sample, we evalu- 
ated the signal to noise (S/N) ratio as the ratio between 
the fluence of the flare and the continuum calculated in 
the time interval [-la",+l<x], where cr is the Gaussian 
width. The minimum detected S/N is 5. Then, to simu- 
late our procedure, we first calculated the median con- 
tinuum light curve from the whole data sample. This 
median light curve at late times is well described by 
a single power law with ^median = IT. On top of 
that we summed a Gaussian flare with the 3 parame- 
ters randomly chosen and uniformly distributed over 
large intervals which fully contain the real data param- 
eter values. From this parent distribution we generated 
a collection of photons. Finally we reconstructed the 
light curve, using the same procedure as we used for 
the real data. In this way we realistically reproduced 
a typical observed light curve. The only significant 
difference is that when we simulated, we assumed a 
continuous observation, whereas the real observation 
is split in different orbits. However, as discussed in the 
previous section, this assumption does not affect our 
conclusions. We repeated the test 14000 times in order 
to have a statistically significant sample of simulated 
light curves. For each randomly generated peak we 
calculated the S/N ratio and we flagged it as identified 
when its S/N ratio exceeded the value of 5 and, at the 
same time, at least 3 points in the light curve lay more 
than 2-cr above the continuum. In Fig. [9] we plot the 
results of the simulations in the (f, At) plane. For each 
(t, At) value we could assign a detection probability; 
the points are the real data. We note that at early times 
our sample data lie in the region with low detection 
probability. This is clearly an effect of the significance 
threshold: at the beginning, the afterglow is brighter, 
the absolute level of the noise is high and a flare can 
be detected only if it is bright enough to have signif- 
icance above the threshold. Given the median contin- 
uum light curve, our simulations show that if flare has 
a ~90% detection probability at 10 ks, at 300s it will 
have a ~ 30% detection probability. At late times the 
simulation results show that the detection probability 
decreases with smaller Af (bottom-right corner of the 
plot): this is an effect of threshold set as the minimum 
number of photons per bin of the light curve. At larger 



times the light curve has a sparse sampling and a faint 
and narrow flare produces only few bins over the con- 
tinuum. Our simulations also show that, in the region 
of the plane defined by At/t > 2xl0~ 3 and by t > 10 4 s, 
the detection probability is uniformly larger than 90%. 
In Fig.|9]we also plotted the line At - t over which we 
do not expect to find any flare. Comparing our sample 
with the simulation probability map we conclude that 
we do not find narrow flares at large t in the areas of 
the parameter plane where we have very high detec- 
tion probability. Therefore, although we cannot evalu- 
ate the completeness of our sample at early times, from 
our simulations we can firmly conclude that the lack of 
narrow flares at late times (typically 10 3 s) is not due 
to incompleteness. 

5. XRT flares vs. BAT pulses 

We investigated if there is a clear link between the 
properties of the pulses detected in the gamma-ray 
burst profile by BAT in the 15-350 keV band and the 
X-ray flares as detected by XRT. 

In order to define a procedure to select and charac- 
terize BAT pulses, we used an adapted ve rsion of the 
criterion defined by Li & Fenimord ( I996I) : we started 
from the 64-ms mask-tagged light curve extracted fol- 
lowing the standard BAT pipeline and searched for 
those bins whose count rates exceed m contiguous 
bins by ncr on both sides. We applied this proce- 
dure with three different combinations of (m, n): (5, 3), 
(3,4), and (1,5) and to all of the curves with mul- 
tiple binning times from 64 ms to 32 s, taking into 
account all of the possible shifts at a given binning 
time. This choice proved to be effective in catching 
different pulses clearly detected by visual inspection. 
We assessed the false positive rate of pulses so de- 
tected with a Monte Carlo test: we took the number 
of 64-ms bins of the longest GRB light curve avail- 
able and simulated 100 synthetic light curves with con- 
stant signal, whose count rates were affected by Gaus- 
sian noise. We applied the same procedure to these 
100 synthetic light curves and found 8 false pulses. 
We then estimated the average false positive rate as 
of 0.08 fake pulses for each GRB light curve. As 
we collected 28 GRBs with a complete BAT light 
curve [GRB 050820A was ignored because Swift en- 
tered the South Atlantic Anomaly (SAA) before the 
gamma-ray prompt emission ceased], we expect about 
2 false pulses. We detected 46 pulses distributed in 28 
gamma-ray profiles, so we can safely assume a negli- 
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gible contamination of the gamma-ray pulses sample 
due to statistical fluctuations. 

Table [4] shows the results of the BAT pulses quest, 
which identified 46 pulses out of 28 GRBs. For each 
pulse, columns 1-6 report: (1) the GRB name it be- 
longs to, (2) the ordinal number of the pulse within 
the GRB, (3) the binning time used to detect the pulse 
(which also corresponds to the uncertainty on the peak 
time), (4) the peak time, (5) the peak rate (counts s _1 ), 
(6) error on the peak rate (counts s _1 ). 

We do not find any clear correlation between the 
number of gamma-ray pulses and the number of X-ray 
flares. Column (1) in Table [5] reports the number of 
gamma-ray pulses found in a given burst; column (2) 
the number of X-ray flares, and column (3) reports the 
number of GRBs with that combination of numbers of 
pulses and flares. The most common case is when the 
burst exhibits one single pulse followed by one, or two 
X-ray flares. 

We tested whether there is any statistical evidence 
that GRBs with many/few pulses are more likely to 
have many/few X-ray flares. Let n y and n x be the num- 
ber of gamma-ray pulses and of X-ray flares of a given 
burst, respectively. We split the sample in two classes 
in two ways: those with n y > 2 ("many pulses"; 
10 GRBs) and those with n y < 2 ("few pulses"; 18 
GRBs); likewise, those with n x > 3 ("many flares"; 1 1 
GRBs) and those with n x < 3 ("few flares"; 17 GRBs). 
From Table [5] one counts 5 bursts with both many 
pulses and many flares. In the assumption of no cor- 
relation between the number of pulses and the number 
of flares, the probability of choosing randomly n > 5 
bursts with many pulses out of 1 1 bursts with many 
flares is about 35%: i.e. given a burst with many flares, 
nothing can be inferred about its number of pulses. 
Similarly, the probability of selecting n > 5 bursts with 
many flares out of 10 bursts with many pulses is 32%: 
i.e. given a burst with many pulses, nothing can be in- 
ferred about its number of flares. We also tried to split 
the sample with different combinations of thresholds 
on n y and n x , but no statistically significant correlation 
has come out. Furthermore, we compared the distri- 
butions of the numbers of pulses derived for the two 
populations, i.e. those with few flares and those with 
many flares. A Kolmogorov-Smirnov (KS) test shows 
no difference between the two subsets, with 88% prob- 
ability that they have been drawn from the same popu- 
lation. Likewise, we compared the distributions of the 
numbers of flares derived from splitting the sample be- 
tween GRBs with few and many pulses, respectively. 



According to the KS test, we cannot reject that the 
two distributions are the same at 99% confidence level. 
We conclude that one cannot infer anything about the 
number of X-ray flares from the number of gamma-ray 
pulses and vice versa. 

We also compared the distribution of the numbers 
of pulses with that of the numbers of flares and a KS 
test does not prove any significantly different origin 
(30% probability of having been drawn from the same 
distribution). 

We also sought any possible correlation between 
the intensity of the pulses and properties of the flares 
as well as between the peak times of either class. To 
this aim, for each GRB in Table [6] we grouped the fol- 
lowing pieces of information: columns 1-3 report the 
GRB name, the number of BAT pulses n y and the num- 
ber of X-ray flares n x , respectively. From column (4) 
up to column (12) the correspondent times are reported 
(referred to the BAT trigger time): the first n y refer to 
the BAT pulses, while the remaining n x refer to the X- 
ray flares. 

For either class we considered those bursts with at 
least two events (i.e., either two pulses or two flares). 
We searched for any correlation between the quies- 
cent time (between two successive pulses, or between 
two flares) and the peak brightness of the following 
event, but our search was unsuccessful. We also stud- 
ied the relation between quiescent time and the ratio of 
the following peak, peak (+1 over the preceding peak, 
peak,. Figure [10] shows two interesting results: firstly, 
there is no clear dependence of this ratio on the qui- 
escent time for both classes. Secondly, the distribu- 
tion of ratios derived from the X-ray flares is consis- 
tent with that of the gamma-ray pulses. In particu- 
lar, if we merge the two sets of ratios, this is consis- 
tent with a log-normal distribution with mean value 
< log (peak ;+I /peak,) >= -0.258 and cr log = 0.68 (see 
Fig. [TTli. If we ignore the two points due to X-ray flares 
with the lowest ratio (see Fig.fTOl. the mean value and 
standard deviation turn out to be -0.157 and 0.41, re- 
spectively (shown in Fig.QJJi. 

We therefore conclude that the relation between 
successive pulses and between successive flares is the 
same: in particular, on average the next event has a 
peak 10 - 0.7 times as high as the preceding, 
while the scatter is between 0.3 and 1.8. This fur- 
ther piece of evidence points to a common origin for 
gamma-ray pulses and X-ray flares. 
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6. Results 

In this section we explore possible correlations be- 
tween the parameters derived in the analysis and sum- 
marize our findings. 

6.1. Gaussian peak time-intensity correlation 

We tested for a correlation between the Gaussian 
peak intensity and the peak position (s since the BAT 
trigger). As shown in Fig.[T2l the correlation is strong, 
with a Spearman rank coefficient r s = -0.539 (num- 
ber of points N = 63, and null hypothesis probability 
nhp= 5.24 x 1CT 6 ). However, it can be argued that this 
correlation is driven by the flares at late times and that 
there is large scatter for t < 10 3 s. In this light, this 
would be an indication that the mechanism producing 
the flares holds no memory of when the trigger time 
occurred. Therefore, the only firm conclusion we can 
draw is that the late flares have a peak intensity which 
is less than the early ones and coupling this with the Af 
results (see 94.21 ) we infer that late flares have a lower 
peak intensity but last much longer so their fluence can 
be very large. 

6.2. EW correlations 

We find a strong correlation between the equiva- 
lent width and the time of the occurrence of the flare, 
fpeak, (r s = 0.729, N = 48, nhp= 4.1 x 10~ 9 ) which 
is mostly due to the large dynamical range in f pea k val- 
ues. Indeed, we find no correlation of EW/t psa ^ with 
fpeak- There is also no correlation between EW/t pe . d k 
and Af /fp ea k (Fig. [T3l which is probably a further in- 
dication that the flares are not related to the under- 
lying continuum and that they originate from the en- 
gine rather than the external shock. We also note 
that EW/ fpeak is generally greater or equal to Af/f pea k 
(solid line in Fig. [T3l ) because the EW calculation is 
sensitivity-limited. The median value of EW/t pe!l k is 
0.5 (mean value 5.7 with standard deviation 25.5). 

6.3. Decay slope-time correlation 

If we consider am as a function of time, we ob- 
tain, for f < 10000 s, that a m = 2.45 + 0.418 f. The 
correlation is only marginal (r s = 0.152, N = 35, 
nhp= 0.382) and a somewhat smaller value is ob- 
tained by using, as stated above, / = 0.05. We con- 
clude therefore that in most cases the exponent of the 
pow er law decay is in agreemen t with the curvature ef- 
fect ( Kumar & Panaitesc ul2000h . In late internal shock 



models, 7q has to be reset every time when the cen- 
tral engine r estarts (IZhang et al. 2006). As shown in 
iLiang et al.l (120061) . if one assumes that the post-flare 
decay index satisfies the curvature effect prediction 
a — /3 + 2, the required To is right before the corre- 
sponding X-ray flares at least for some flares. This 
lends support to the curvature effect interpretation and 
the internal origin of the flares. In a few flares, how- 
ever, the giant flare observed in GRB 050502B being 
the best example, the deca y slope is much steeper if 
To is put near the peak (see, Dermerll2004 ; Liang et al. 
l2006h . 

6.4. Afftu/Afrise — T90 correlation 



Du ring the prompt emission, as tested by lNorris et al 
(11996b . shorter bursts tend to be more symmetric and 
the width of the burst tends to be correlated with 
Aff a u / Af r i S e in the sense that longer bursts tend to 
show a larger ratio, or Aff a u/Af r i se ~ 2-3, a value 
that agrees quite well with the mean (A?faii/AZri se ) = 
2.3 5. This effect has bee n quite clearly simulated 
by iDaigne & Mochkovitchl dl998l) . The flare sam- 
ple was used to test for this effect. We used T90 as 
a reference time to minimize the bias we may have 
in the curve subtraction when the signal of the flare 
is weak. In addition, we considered both expressions 
such as F(t) = a(\ - bx~ c ) and simple power laws 
(F{t) = k x~ m ) to model the sides of the flares. Using 
/ = 0.05 the difference in the width (T90) of the flare 
inferred from the two fits is negligible for the scope of 
this work. 

As shown in Fig. [14] we find a tentative corre- 
lation between the ratio Aff a n / Af,-i se and T90 (r s = 
0.543, = 24, null hypothesis probability, nhp= 
6.15 x 10~ 3 ). Such a correlat ion was pointed out by 
Daigne & Mochkovitchl (119981) in their simulations of 
the prompt emission. 

6.5. Summary of Results 

We gathered a sample of 33 light curves drawn from 
all GRBs detected by Swift, INTEGRAL and HETE-2, 
which had an XRT follow-up and which showed either 
large-scale flaring or small scale (mini-flaring) flick- 
ering activity. None of the INTEGRAL- or HETE 2- 
triggered bursts showed any flares (however, note that 
these burst were observed by XRT much later than the 
Swift-triggered ones). For 30 of these bursts, we per- 
formed a full statistical analysis, by fitting the contin- 
uum light curve beneath the flares (the XRT canonical 
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light curve shape) with a multiply-broken power law 
and the flares with a sample of analytical functions. 
Our sample of Gaussian fits consists of 69 flares, for 
48 of which we calculated EWs by numerical integra- 
tion, for 35 we could determine a decay slope, and for 
24 of them T90, Aff a u/A?ri se and Af/f. Our results can 
be summarized as follows. 

1. Flares come in all sizes and shapes and can 
be modelled with Gaussians (symmetrically 
shaped) superimposed on a multiply-broken 
power-law underlying continuum. However, for 
a more accurate description, in many instances 
an exponential rise followed by a power law de- 
cay or power law rise followed by a power law 
decay is required to produce good fits. 

2. Flares are observed in all kinds of GRBs: long 
(32 GRBs) and short (2 GRBs), high-energy- 
peaked or XRFs (32 vs. 2); they are found both 
in early and in late XRT light curves. 

3. The equivalent widths of our sample, which 
measure the flare fiuence in terms of the under- 
lying continuum, range between 8 s and 7x 10 s s. 

4. The distribution of the ratio Af/f, as defined by 
the width and peak of the Gaussians flare mod- 
els, yields (Af/f) = 0.13+0.10. Our simulations 
show that our time resolution allows us to sam- 
ple flares that may have Af/f < 0.1, so that the 
above values are not the result of the biases in 
our sample or our fitting procedures. Our simu- 
lations also show that there are no sharp (small 
Af/f) flares at large times. 

5. The decay slopes ar&u range between 1.3 and 6.8 
and generally agree with the curvature effect. 

6. The ratio of decay and rise times range between 
0.5 and 8. 

7. Correlations are found between 

(a) fp ea k-peak intensity (strong); 

(b) EW-fpeak (very strong); 

(c) ff fa ii-f (poor); 

(d) AWA?ri Se -T 9 o (tentative). 

8. We do not find any clear correlation between the 
number of gamma-ray pulses and the number of 
X-ray flares. One cannot infer anything about 



the number of X-ray flares from the number of 
gamma-ray pulses and vice versa. We also con- 
clude that the relation between successive pulses 
and between successive flares is the same: in 
particular, on average the next event has a peak 
jq-0.157 ^ o.7 times as high as the preceding, 
while the scatter is between 0.3 and 1.8. This is 
a piece of evidence pointing to a common origin 
for gamma-ray pulses and X-ray flares. 



7. Discussion 



The analysis of the flares in the present sample 



light curve (|Chincari 


tii et al.ll2005 


; Nousek et al.l2006: 


O'Brien et all 2006; 


Zhane et alj 


2006) make it clear 



that the onset of the XRT observation corresponds to 
the late tail of the prompt emission as defined in the 
current model. Flares are often observed in the early 
XRT light curves. Their slopes do not conflict with 
the curvature effect limit; they simply nee d a different 
interp retation and a proper location of Tq dLiang et al. 
l2006h . 

A similar reasoning explains the decay slope of the 
flares. We have se en, in agreement with the finding of 
Liang et al. (2006), that the decay slope is very sensi- 
tive to the definition of To and that if this is located at 
the beginning of the flares we are within the constraint 
of the curvature effect. This essentially means that the 
shock, after reaching the maximum luminosity, is not 
fed anymore and fades out. Some of the uncertain 
or critical cases of flares may be due to the presence 
of blends. Blends and superimposed mini flares are 
indeed very common and we can observe them very 
clearly in all those cases in which the statistics are very 
good. Although the analysis may be affected in part by 
this contamination, the results remain robust. Indeed, 
the contamination makes our results even more robust 
since the detection of unseen blends would make the 
selected Tq large, thus decreasing the measured slope 
and width of the flares. 

We also considered the possibility of a correlation 
between the characteristics of the prompt emission as 
observed by BAT and the frequency of flares detected 
by XRT. We found no correlation. This simply means 
that the flares are random events and are not related to 
the way the prompt emission develops in time. For in- 
stance, there could be an initial flickering, due to the 
collision of highly relativistic shells followed by ran- 
dom flare events due to the collision of slower resid- 
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ual pellets, as discussed below. The contamination to 
our sample due to the fact that some of the early XRT 
flares are the tail of the late prompt emission does not 
change this result. However, this needs to be further 
investigated using a larger statistical sample. 

Furthermore, we have shown that our analysis is not 
affected by bias in the detection of high-intensity late 
flares and that such flares never show a peak of inten- 
sity as strong as those observed in the early flares. On 
the other hand, due to their rather long duration, these 
flares are also very energetic. 

Most of the indications we have so far seem to lead 
toward an activity that is very similar to that of the 
prompt emission, with flares that are superimposed on 
a very standard light curve. This has been observed 
both in long and short bursts. 



In light of the calculations of lloka et al.1 (2005), 



we calculated AF/F and At/t values from our flare 
sample and plotted them over the kinematic ally al- 
lowed re gions for af t erglow variabilities, as shown in 



Ioka et al.l (12005b distinguish between four 



Fig. El 

cases: (a) dips, arising from non-uniformity on the 
emitting s urface indu c ed, e.g ., by density fluctuations, 
[eq. (4) in lloka et al, (2005)] ; (b) bumps due to den 



sity fluctuations dWang & L oeb 2000; Lazzati et al 



20021 lDai&Lull2002h [eq. (7) in lloka et all (120 05)1: 
(c) bumps due to patch y shells (uvleszaros et al.lll998 : 
Kumar & Piranll2000ah . for which At > t; (d) bumps 
due to refreshed sho c ks dRees & Meszarosl f l998: 



Panaitescu et all [1998; Ku mar & Piran 2000b), for 



which At > t/4. 

Our finding s are c onsist ent with the co n clusio n 
of IZhang et all d2006l) and lLazzati & Pernal d2007h . 
the latter based on a prelim inary presentation of our 
dataset in Chincarini 00061) . i.e., a sizable fraction of 
the flares cannot be related to external shock mech- 
anisms. In particular, only one point (correspond- 
ing to a flare in GRB 051117A) lies in the region of 
At > t, where flares are consistent with the patchy 
shells model; only three points (including the the early 
flare of GRB 050502B) lie in the region of flares that 
can be caused by ambient density fluctuations; finally, 
only 29/69 lie in the region that describes flares due to 
refreshed shocks. Among the rest, 10/69 can only be 
due to internal shocks. 

Pern a et al. I d2006l) proposed that X-ray flares are 
due to accretion of a fragmented disk. Due to vis- 
cous evolution, blobs far from the central black hole 
takes longer to be accreted and are therefore more 



spread-out when accretion occurs. The accretion rate 
is correspondingly lower. This naturally gives a peak 
luminosity-flare epoch anti-correlation as has been re- 
vealed by the data. This same merit could be re- 
tained if a magnetic barrier modulate a continuous ac- 
cretion flow near the black hole at different epochs 
dProga & Zhangll2006l) . 
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Table 1 

grb xrt light curve sample. 
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w > 




(1) 


050406 




5 ± 1 


9.0 x 10~ 8 




1 


XRF 


050421 




10.3 + 2 


(1.8 +0.7)x 10~ 7 




2 




050502B 




17.5+0.2 


(8.0+ 1.0) x 10~ 7 




3 




050607 




26.5 


(8.9 + 1.2) x 10~ 7 




4 




050712 




48 ±2 


1.8 x 10~ 6 




5 




050713A 




70+10 


(9.1 +0.6)x 10~ 6 




6 




0507 14B 




55. 


(6.5 + 1.4)xl0- 7 




7 


XRF 


050716 




69 ±1 


(8.3 + 1.3)xlO- 6 




8 




050724 


0.258 


3 + 1 


(6.3 + 1.0) x 10- 7 


9 


10 


Short 


050726 




30. 


(4.3 + 0.7) x 10~ 6 




11 




050730 


3.967 


155+20 


(4.4 + 0.4) x 10~ 6 


12 


13 




050801 




20 ±3 


(4.4+ 1.0) x 10~ 7 




14 




050802 




13 ±2 


(2.8 +0.3)x 10~ 6 




15 




050803 


0.422 


85 + 10 


(3.9 + 0.3) x 10- 6 


16 


17 




050814 




65 +4 ° 

—20 


(2.17 + 0.36) x 10- 6 




18 




050819 




36 + 4 


(4.2 + 0.8) x 10~ 7 




19 




050820A 


2.612 


26 ±2 


(1.9 + 0.2) x 10~ 6 


20 


21 




050822 




102 + 2 


(3.4 + 0.3) x 10- 6 




22 




050826 c 




35 ±8 


(4.3 + 0.7) x 10~ 7 




23 




050904 


6.29 


225 ± 10 


(5.4 + 0.2) x 10~ 6 


24 


25 




050908 


3.3437 


20 ±2 


(5.1 +0.5) x 10~ 7 


26 


27 




0509 15A 




53 ±3 


(8.8+0.9)xlO- 7 




28 




050916 




90+10 


(1.1 + 0.4)x 10~ 6 




29 




050922B 




80+10 


(1.8 +0.3)x 10~ 6 




30 




051016B d 


0.936 


4.0 + 0.1 


(1.7 +0.2)x 10~ 7 


31 


32 




051117A 




140 + 10 


(4.6 + 0.16) x 10" 6 




33 




051210 




1.4 + 0.2 


(8.3 + 1.4)x 10~ 8 




34 


Short 


051227 




8.0 + 0.2 


(2.3 +0.3)x 10~ 7 




35 




060108 




14.4+1 


(3.7 + 0.4) x 10~ 7 




36 




060109 e 




116 + 3 


(6.4+ 1.0) x 10~ 7 




37 




0601 11A 




13 + 1 


(1.18 + 0.05) x 10- 6 




38 




060115 


3.53 


142 + 5 


(1.8+0.2)xl0- 6 


39 


40 




060124 f 


2.296 


321+2 


(1.40 + 0.03) x 10- 5 


41 


42 





a GRBs with number in italic were considered for their behaviour, but did not offer sufficiently 
high statistics to allow full analysis (see 

b Drawn form refined BAT GCN Circulars in the 15-150keV band. 

C A low-signal late-time flare is observed and no analysis was performed. 
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d A flattening in the XRT light curve is observed starting from t ~ 200s and lasting through the 
first SAA data gap. A fit with a Gaussian centered at t ~ 650s provides a significantly worse fit than 
a combination of power laws, hence this event was not included in the restricted sample. 

e A flattening in the XRT light curve is observed starting from t ~ 10 3 s and lasting through the 
first SAA data gap. 



f As reported in iRomano et alj d2006al) . a separate fit was performed to the prompt and the after- 
glow parts of the X-ray light curve. Here we do not consider the spikes in the prompt. 



Refer ences. — (1) [ Rrimm et aL ,. 
d2005bh: (4) iRetter et al| d2005h: (5) 



(2005a); 
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(2005 bl): (2) ISakamoto et al.1 d2005al): (3) ICummings et al 

Palmer et al.1 d2005ch: (7) 



. , Markwardtet al , 
Tueller et all d2005al): (8) IBarthelmv et al.1 d2005bl): (9) IProchaska et al.1 d2005bl): (10) iKrimm et al.. 
d2005ah: (ll)lBarthelmy et al.1 d2005dl): ( 12) l(Ihen et alJ d2005b: (13)lMarkwardtetalJd2005cl): (14) 



Sakamoto et alj d2005ch: ( 1 5) IPalmer et al] d2005al); (16) iBloom et all d2005l): (17) IParsons et al 
d2005h: (18)lTueller et al.l d2005bl) : (19) IBarthelmv et all d2005ch: (20) IProchaska et all d2 005a): 



(21) ICummings et all d2005ah: (22) iHullinger et al I d2005ah: ~C23) iMarkwardt et all d2005bh: (24) 
iaislip et all d2006b: (25) ISakamoto e t alJ d2005bl): (26) iFueazza et al.1 d2005l): (27) ISato et al ' 



(2006): (42)|Romano et al.l d2006al) . 



(2005aT); (35)lHumnger et all d2005bl): (36) ISakamoto "et al] 120061) : (37) Palmer et all d2006h: (38) 
Sato et al. (2006); (39) Pirano monte et al.l OOOd) : (40) lBarbier et al.ld2006h : (41) lMirabal & Halpernl 



2005bl): (28) IBarthelmv "et al.l d2005al): (29) iFenimore et al.l d2005b: (30) Hullinger et all d2005 
(3 1) ISoderberg et all d2005l): ( 32) iBarbier et al.l d2005b: (33) IPalmer et alj d2005bl): (34) ISato et al 



14 



Table 2 

Fits to the XRT light curves: continuum parameters. 



GRB 




?bl 


«2 a 


tb2 


a 3 d 






(s) 




(s) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 



050406 

050421 

050502B 

050607 

050712 

050713A 

0507 14B 

050716 

050724 

050726 

050730 

050801 

050802 

050803 

050814 

050819 

050820A 

050822 

050904 

050908 

0509 15 A 

050916 

050922B 

051117A 

051210 

051227 

060108 

0601 11A 

060115 

060124 c 



1 - 58 -ai7 
3 10+ 011 

J ' 1 -0.09 

75+ 004 

u - -0.04 

1 65+ 017 

-Q16 

2 17+°^ 

-0.67 

7 - i6 -S 

6 79 + " : ° 5 

u - -0.38 

1 09+0.02 

-Q.Q7 
1 co+0.07 
• J -Q.Q7 
+0.04 

8+ 



m 

0.16 



0.95 
0.28 
0.57 

0. 27 b ' 

4 54+0.26 

-0.29 

3 26+°- 12 

-0.20 

3.22 b 
9 9C+0.14 

-0.17 

2.99 b 

1 57+°- 12 

1. J/_g 13 

j j 2+0.06 

42+^ 

u ' HZ, -0.27 

95+ 030 

T OT+0.37 

066+°^ 
-Q10 
9 £ro+0.25 

2 50+ 015 
2 6O+ 055 

_4 9C+O.36 
+•^-0.44 
o 9Q+O.2I 

-0.29 

44+0.07 

-0.08 



(4.36^ |) x 10 3 

(15.2+5-2) x 10 4 
1.45 x 10 3b 
3.44 x 10 2b 
1.12x 10 2b 

(3.90+ / 31 ) x 10 2 

(4 - 70 -Q15 ) x 102 
(1.90 

(8.53+| 32 ) x 10 
(5.52+of 4 ) x 10 3 



+0.05% x 1 r»2 
+I.32\ w in3 



(4.67+f|) x 10 2 
(8.80+f '°) x 10 3 
(4.46+||) x 10 2 
(9.99_« : ») x 10 2 

(8.18 x 10 2b 
(2.00+°j 4 ) x 10 2 

7.50 x 10 2b 
(3.35^) x 10 2 

(1.74+2-36) x 1Q 3 



7.37 x 10 2b 
(2.54+»^) x 10 2 
(3.25+°f>) x 10 2 
(5.74+$|) x 10 2 
(1.0-1 1.5) x 10 3 



50+ 013 



1+0.14 
-0.25 



0.52 
3.12 
0.81 b ' 



49+0.10 
U - H -0.09 

O + 1.40 
8 - 8 -1.6Q 
c O+0.50 
J -0.30 

1 89+°- 16 

1,0 -0,20 
J 97+O.O6 

1 94^:8^ 

1 -0.08 

1.60+ 019 

03+^ : ^ 
56+ 009 

U - J "_0.14 

0.27 b 
0.03 b 
0.40 b 

2 26+ 011 
•^"-o.io 



1 20+ 020 
1 - z,u -o.io 



0.18 b 

37+0- 05 

6 2**® 

"•^"+0.28 

70+ 005 

u - -0.08 

1.21 + 0.04 



1.54 x 10 4b 
8.39 x 10 2b 

(8.03^1') x 10 



(5.53_°j°) x 10 2 



(1.27!« |) x 10 4 
(8.46!j; 24 ) x 10 4 



3.10x 10 3b 
(1.87+°- 36 -) x 10 3 
(7.38+jjfj) x 10 2 
(3.9 1+ 2 ] 3 ) x 10 4 
(1.05!° : ! 7 )xl0 5 



j 34+0.39 
043+8$ 

"•^ -0.29 

79+ - 33 

u - -0.32 



0-78!!J;! 3 



1 cq+o.03 

,J -0.04 
9 44+O.34 

^■•++-0.46 



(4.79^ 2 )xl0 3 1.27+oos 

(2.22xl0 4b 1.72 b 
(1.70^; 49 )x 10 4 0.50 b 



1.22 b 

1 99+O.O9 
-0,09 

90+ 05 

1 o 1 +0.22 
1,J1 -0.20 

1.58 + 0.06 



a These slopes do not strictly correspond to phases I, II, and III of the canonical XRT light 
curve. 

b Parameter fixed. 

The fits of prompt (first orbit) and afterglow were performed separately. The first break (fbi) 
is not well defined, since it occurs during a SAA passage, that lasts from ~ 1000 to ~ 1 1500 s. 
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Table 3 

Fits to the XRT flares: Gaussians, power-laws and exponentials. 



GRB 




Gaussian 




AF/F 


EW 


am 


T90 


Alfall/A/nse 


At// 11 




Center (s) 


Width (s) 


Norm (count s 1 ) 




(s) 




(s) 






(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 



050406 
050421 



050502B 



050607 
050712 



050713A 

050714B 
050716 

050724 

050726 
050730 



050801 
050802 
050803 



050814 
050819 
050820A 
050822 



050904 



050908 



050915A 
050916 



050922B 



051117A 



211+5 

m+° 

154+| 
719+1 
3343 1 +4057 

74637+^' 

330+ 8 
245.7+" 
486. 1+™ 
913-5™ 
112.2+"° 
173.4+]]! 
399.S+?]? 
126.2+3] 
399+*' 
175+° 
382+= 
275*1 
327™ 
(5.7+°]) x 10 4 (1 
168+= 
273+J. 
131.8+>" 
234.2+IJ 
436.5+H 
6S5.8+ 2 -' 8 , 
742 ' 
4526.2+!^;* 
10223 6 +2l)! 
12182.9 
284+* 
464+ 
332+™ 
604 
1201 
2286+™ 
177+ 7 

1 ' '-19 



241+ 
142 - 7 -u 

9/11 S+W 



241. i 

465 - 7 -i. f , 
448.6+3' 
975.5+g| 
1265.5+1*° 
71 12.832*3 
16682.2+1?* 
31481.1+™ 
146!'" 
425+jf 
107+| 
18750+™' 
21463+™' 
375+5 
490+ 8 
858+]° 
132+5 



17.9+] 2 , 3 

1 7+°' 
'■ '-0.1 

6.2+« 

ioo.i!|j 2 

6273.0+»| 
26742.9!™;' 
36.1+g 
31.1+|g 
16.7!+' 
98 

-0.5 

16.3+?;° 

23 3+ 5 -' 
-3-S 
10.8!] | 

52.4!" 

48!] 

16.3!; 

30.6!"' 

12 7+" 
lz - -5.0 

.9+°])x 10 4 
8.2!$] 
27.0!=] 

^>-> 7+24.4 

-5,3 

14.5+3 
38.5+lf 
23.8+]]] 
10 

408. l+_lff 
847+ 3 
383.2 

49 <+42.5 

-42,5 

29.0!g 
189.2 



299 n+ 423 - 7 

13.9+]' 5 72 

9 5+03 

15 2+'° 

1J - -0.8 

12 4+'-' 

49 0+ 2 - 3 

45 9+ 4 - 5 

62.8+*]? 

81.6+]° : 2 

790 4 +1 ° 3J 
-!i,5 
3194 8+ 2122 

3 "+-°-227.i 

7149 6+<>87.4 

' '+'■ -647.6 

23 +-23 

45+' 8 
15 5+ 5 ,s 
425+205 

2222+® 2 

zzzz -360 

9 2+ 2 -' 
37.7+1 

123 8 
48!] 



4 6+ 12 

191 9+1186.8 
-107 .8 

4.7+f 4 



012+0006 

0^027+IP 
15.8!];] 



7.3 
7.4 



+1:5 



170.5+$ 
23 5+ i4 
24.9+ 3 - 4 , 



55.6+, 



°-9.6 
4 4+0.8 
+•+-0.9 

9.0!] 4 

3-8-1.3 
7.2!!] 
3.1+to 

030 +oJ ° 04 

J1 -1.6 

6.6!] 

' -0.8 
9.0+0.6 

5 ig+0'69 

3.0 

°- 87 -„:!„ 

0.4 



2 14+0.46 

0.8+ff 
1.00 
0.67 
12+ 002 

2 ,+-°.o° 2 

231.0+° 



54.7 



+f 
15.5+1 
43 5+138 

+ J - J -1.46 

20.7!!] 
1 +oii 

'•"-0.2 

1 l+°- 4 
'■'-Q.3 

1 6 +0J 
n 77+0.05 

"•"-0.04 

0.31!°°; 
7 17 +8:^3 

24+P:? 7 

12.20!]] 

0.2+°,!' 

o.i!™ 

23.0!] 
6.7+ 11 
22.0" 
102+5" 



7.71 


686 




184.0 


1.520 


0.882 


12.35 












0.84 












38.55 


127320 


6.3 ±0.38 


523.6 


1.450 


1.352 


0.28 












1.57 


432630 


4.67 ± 0.34 








21.84 


1813 


3.39 ± 0.24 


266.5 


1.610 


0.798 


1.37 




... a 








5.10 


165 


2.87 ± 0.41 








6.35 












19.77 


190 


2.92 ± 0.25 


49.5 


3.070 


0.445 


3.97 


94 


3.1 ±0.2 


82.9 


2.790 


0.494 


8.27 












7.26 












69.86 


2313 


3.2 ± 0.43 


344.8 


3.410 


0.928 


0.48 


240 


0.51 ±0.24 


622.0 


4.750 


3.514 


0.57 


383 


2.13 ±0.51 


482.9 


2.700 


1.283 


1.35 


84 


2.52 ±0.5 








1.58 


67 


4.43 ±0.8 








11.99 


737109 


3.13 ± 0.19 


112365 


1.720 


2.045 


0.46 


8 


3.7 ±1.2 


33.0 


0.492 


0.199 


1.57 


126 


3.5 ±0.53 


122.0 


1.120 


0.446 


1.33 












1.00 


43 


4.9 ±1.1 








2.02 


370 










1.32 


224 










0.78 












0.37 


350 










1.34 


1897 










0.87 












0.91 












2.25 


159 


2.54 ± 0.35 


926.3 


5.300 


2.327 


0.85 


65 


1.7 ±0.8 








4.05 


357 


3.1 ± 1.2 








2.66 


404 


1.3 ±3.2 








1.30 












0.67 


8 










77.45 












1.09 


59 


4.34 ±0.17 








1.50 


129 


2.78 ±0.17 


110.3 


1.280 


0.459 


29.89 


6851 


5.06 ±0.18 


328.0 


0.630 


0.708 


2.22 


401 


4.52 ± 0.32 








0.62 


162 










1.23 


364 










88.46 












292.29 












166.67 












1.72 


88 


b 








6.29 


1132 


2.36 ±0.11 


295.6 


2.660 


0.727 


13.36 


43 


3.35 ± 0.38 








25.22 


130717 










14.34 












0.97 


221 


1.66 ±0.33 


175.4 


1.330 


0.466 


0.69 


410 




254.7 


2.020 


0.508 


14.64 


14336 


6.76 ± 0.42 


464.9 


1.420 


0.572 


3.23 


27 


2.72 ± 0.46 


331.9 


2.5 


2.192 
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Table 3 — Continued 



GRB 


Center (s) 


Gaussian 
Width (s) 


Norm (count s 1 ) 


AF/F 


EW 

(s) 


<*fall 


T 90 
(8) 


Af fall /A'me 


At/t* 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 




376^ 


203+" 


4T_l 


3.01 


195 












955^, 


69+« 


29+> 


3.41 


395 


3.51 ±0.28 










1110+5 


50+J 


27+2 


3.51 


1201 












1341+5 


43+3 
135+] 


49+3 


7.21 






603.4 


8.060 


0.453 




1516+^ 


30+6 


4.78 




6.6 ±4.0 








051210 


134.435 
216.2 


10.4+1] 
63.1 


4 7+1.7 
^■'-1.7 

0.62 


1.18 
0.53 


16 


4.05 ± 0.57 


49.2 


0.490 


0.360 


051227 


124.2 


10.5 


5.15 


0.88 




2.05 ± 0.5 








060108 


303.5+|« 


44 5+125.5 
++ J -30.5 


0-3+g;!I 


1.83 


25 










0601 11A 


95.1+1* 


22.8+fi 
18.4^1 


67.6+IJ 


165.51 


73 


3.53 ±0.39 


144.4 


0.800 


1.405 




166.9+'* 


34.7+IJ 


7.78 


54 


4.5 ± 1.2 


120.7 


1.230 


0.719 




280.1+H 


20 6+ 173 

/:? - 1 -23.5 


2.11 


931 


6.51 ±0.4 


177.5 


2.430 


0.620 


060115 


431.9 + [|| 


1 gi+o:5J 


2.53 


144.2 











''Using At = (the time defined in terms of f = 0.05, j |4.3l and t = ?p ea k- 

b GRB 050712 and GRB 050908 have a first flare that quite likely is part of the prompt emission. In addition the decay does not show a 
very high statistics. 
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Table 4 

Properties of the 46 gamma-ray pulses detected from the BAT light curves of 28 GRBs with X-ray flares. 



GRB 


N 


Bin 


Peak 


Peak 


Error on Peak 


Name 




T(s) 


Time (s) 


Rate (counts s _1 ) 


Rate (counts s~') 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


050406 


1 


O 1 "7 

2.176 


O O A 

2.24 


A A A 1 A A 

0.04144 


A AAC A £ 

0.00546 


AC A A O 1 

050421 


1 


10.688 


1 1 ATI 

1 1.072 


A A 1 COO 

0.01538 


A AAOC/C 

0.00256 


AC AC AO D 

050502B 


1 


A 1 O /I 

0.384 


A OC A 

0.864 


0.23168 


A A 1 C O 1 

0.01581 


05060 / 


1 
I 


1.216 


1.656 


A 1 OAO£ 

0.12026 


A AAA 1 A 

0.00919 




2 


13.952 


1 c ooo 

16.888 


A AOTO 1 

0.03/21 


A AAOC/C 

0.00256 


AC AT 1 O 

050712 


1 


OA A £. A 

30.464 


O/C A 1 O 

26.912 


A AT AAC 

0.03095 


A AAOTT 

0.00237 


AC AT 1 1 A 

0507 13A 


1 


C 1 1 O 

5.312 


C /I OA 

-54.89 


A A A O A A 

0.04249 


A AA/CTC 

0.00675 




2 


1.664 


O T 1 O 

2.712 


0.53963 


A A 1 TO A 

0.01724 




3 


£ A AA 

6.400 


1 A O /I 

10.84 


A /I AC/: 1 

0.40561 


A AAAO/C 

0.00926 




A 

4 


1 CO A 

3.520 


£A /I 

69.4 


A ACCOC 

0.05585 


A AAC A O 

0.00542 




c 

J 


1 o r\oo 
12.052 


1 1 /T TA 
1 16. /0 


A A1 OA A 

0.01849 


0.00233 


AC AT 1 A D 
050/ 1413 


1 

1 


oo /; o 
23.68 


CO OAO 

52.392 


A AO£ 1 T 

0.0261 / 


0.00313 


AC AT 1 /; 

050/ lo 


1 

1 


/I TOO 

4.288 


11 O 1 

11.81 


A OOO 1 A 

0.22210 


A A1 1 CO 

0.01 168 




2 


1 1.648 


46.94 


A lOTCO 

0.12358 


A AAjJ/I O 

0.00648 


050/24 


1 


A 1 OO 

0.128 


A 1 A/1 

0.104 


1 OAA/1T 

1.2994/ 


A ATA/; A 

0.0/964 




2 


A A£ A 

0.064 


O 1 A AO 

210.92 


A OAAOT 

0.2993 / 


A A/1 OO 1 

0.04221 


050726 


1 


2.368 


1 to ot 

-173.87 


A AAC1/C 

0.09576 


A A 1 C AC 

0.01595 




2 


i~i 1 A/I 

7.104 


T OA 

7.89 


A f\A A A 

0.10490 


A AAA/; O 

0.00968 


AC AT OA 

050/30 


1 


1 O OO /I 

13.824 


n /too 
1 /.432 


A A/1 cno 

0.046/8 


A AAO O C 

0.00286 


AC AO A 1 

050801 


1 


A C 1 O 

0.512 


A CAO 

0.592 


A OAOOO 

0.20882 


A A 1 A TO 

0.014/2 


AC AO AO 

050802 


1 


1 O TJ /T 

12.736 


1 O AA 

13.90 


A AO C A A 

0.02504 


A aaoc/c 

0.00256 


AC AO AT 

050803 


1 


1 (CAA 

1.600 


1 1 /TO 

1.168 


A OO A £0 

0.32463 


A AOO CT 

0.02357 


ACAO 1 A 

050814 


1 


16.512 


18.856 


A A /I /C/C 1 

0.04661 


A A A A C A 

0.00450 


ACAO 1 A 

050819 


1 


1 no 

16.128 


23.224 


A AOO /I 

0.02246 


A AAOO 1 

0.00231 


050822 


1 


3.264 


3.912 


A 1 CT 1 A 

0.15 / 19 


0.012/2 




2 


0.96 


48.52 


0.242 16 


0.01586 




3 


4.352 


60.04 


0.10462 


0.00577 




4 


2.624 


103.56 


0.04600 


0.00684 


050904 


1 


6.976 


29.768 


0.06191 


0.00493 




2 


15.808 


125.128 


0.05769 


0.00229 


050908 


1 


3.072 


3.776 


0.07615 


0.00586 


0509 15 A 


1 


5.376 


5.496 


0.05597 


0.00436 




2 


0.768 


14.584 


0.11584 


0.01225 




3 


1.920 


44.6 


0.03813 


0.00520 


050916 


1 


16.448 


51.304 


0.03522 


0.00376 


050922B 


1 


15.168 


52.072 


0.06915 


0.00699 




2 


14.336 


103.4 


0.05754 


0.00705 




3 


1.408 


263.464 


0.06536 


0.00757 




4 


1.216 


271.656 


0.06913 


0.00825 



18 



Table 4 — Continued 



GRB 


N 


Bin 


Peak 


Peak 


Error on Peak 


Name 




T(s) 


Time (s) 


Rate (counts s _1 ) 


Rate (counts s _1 ) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


051117A 




8.704 


11.264 


0.08273 


0.00388 


051210 




0.640 


0.88 


0.1057 


0.01215 


051227 




0.640 


0.80 


0.12975 


0.01213 


060108 




2.56 


3.304 


0.08358 


0.00554 


0601 11A 




2.816 


5.792 


0.19368 


0.00639 


060115 




5.312 


6.352 


0.05227 


0.00438 




2 


3.52 


98.192 


0.09428 


0.00513 
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Table 5 

Frequency distribution of the number of BAT pulses vs. the number of X-ray flares in 28 bursts. 



Number of BAT Pulses Numer of X-ray flares Frequency 



(1) 


(2) 


(3) 




1 


8 




2 


4 




3 


4 




6 


1 




8 


1 


2 


1 


2 


2 


2 


2 


2 


3 


1 


2 


6 


1 


3 


1 


1 


4 


3 


2 


5 


4 


1 
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Table 6 

Central times of the y-RAY pulses and X-ray flares for each GRB. 



GRB 


BAT 


XRT 










Times 










Name 


«7 


tlx 


(s) 


(s) 


(s) 


(s) 


(s) 


(s) 


(s) 


(s) 


(s) 


to 


(2) 


(?) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


050406 


1 


1 


2.24 


211.0 
















050421 


1 


2 


11.072 


111.0 


154.0 














050502B 


1 


3 


0.864 


719.0 


33431. 


74637. 












050607 


2 


1 


1.656 


16.888 


330.0 














050712 


1 


3 


26.912 


245.7 


486.1 


913.5 












050713A 


5 


4 


-54.89 


2.712 


10.84 


69.4 


1 16.70 (D) 


112.2"" 


126.2 W 


173.4 


399.8 


0507 14B 


1 


1 


52.392 


399.0 
















050716 


2 


2 


11.81 


46.94 


175.0 


382.0 












050724 


2 


3 


0.104 


210.92 


275.0 


327.0 


57000.0 










050726 


2 


2 


-173.87 


7.89 


168.0 


273.0 












050730 


1 


8 


17.432 


131. 


234.2 


436.5 


685.8 


742.0 


4526.2 


10223.6 


12182.9 


050801 


1 


1 


0.592 


284.0 
















050802 


1 


1 


13.90 


464.0 
















050803 


1 


3 


1.168 


332.0 


604.0 


1201.0 












050814 


1 


1 


18.856 


2286.0 
















050819 


1 


1 


23.224 


177.0 
















050822 


4 


3 


3.912 


48.52 


60.04 


103.56 


142.7 


241.8 


465.7 






050904 


2 


6 


29.768 


125.128 


448.6 


975.5 


1265.5 


7112.8 


16682.2 


31481.1 




050908 


1 


2 


3.776 


146.0 


425.0 














0509 15A 


3 


1 


5.496 


14.584 


44.6 


107.0 












050916 


1 


2 


51.304 


18750.0 


21463.0 














050922B 


4 


3 


52.072 


103.4 


263.464 


271.656 


375.0 


490.0 


858.0 






051117A 


1 


6 


11.264 


131.7 


375.9 


955.0 


1110.0 


1341.0 


1516.0 






051210 


1 


2 


0.88 


134.4 


216.2 














051227 


1 


1 


0.80 


124.2 
















060108 


1 


1 


3.304 


303.5 
















0601 11A 


1 


3 


5.792 


95.1 


166.9 


280.1 












0601 15 


2 


1 


6.352 


98.192 


431.9 















'The BAT pulse occurs simultaneously with the two X-ray flares. 
'Errors on columns 4-12 are the binning times (see, $5\. 
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Time since BAT trigger (s) Time since BAT trigger (s) 



Fig. 1. — Flare fits. The thick line is the best fit to the XRT data (filled circles) with a (multiply-)broken power law 
plus a number of Gaussians (see Table|3]for the fit parameters). The continumm and Gaussian parameters are reported 
in Tableland Table [3] respectively. For GRB 060124 we considered the prompt and afterglow portion of the light 
curve separately. 
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Fig. 1. — Continued. 
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Fig. 1. — Continued. 
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Fig. 1. — Continued. 
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1 2 3 4 5 

log t peak 



Fig. 3. — Distribution of the peak times of the flares in excess of the canonical XRT light curve. The times are referred 
to the trigger time, and are not corrected for redshift. 
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-3 -2-10 1 

log (At/t) 

Fig. 4. — Distribution of the ratio of the flare duration vs. the time of occurrence At/t, obtained fitting the flares 
with Gaussian models ( jH4.2l i. where Af is the width of the Gaussian and t is the Gaussian peak time. This ratio is 
independent of redshift. 



28 



GRB050502B 



a- 4 




X 



















































































» 




1 

• 

1 










1 



200 400 600 800 1000 1200 1400 
Time Since Trigger 



GRB060111A 




250 300 350 400 450 500 
Time Since Trigger 



GRB050730 




350 400 450 500 550 600 650 
Time Since Trigger 



GRB050730 





•A 
» } 


i 








1 


I • 


n 










\i 
t\ 










T 


it 


• 

• 













4000 6000 8000 10000 12000 
Time Since Trigger 



Fig. 5.— Different flare morphologies as represented by GRB 050502B, GRB 050730, and GRB 0601 11 A. For 
GRB 050730 different flares are best fit by different laws (two power laws for the first, and an exponential rise followed 
by a power-law decay for the second one.) 



29 



8 




30 



2 4 6 8 10 

At-fall AW 

Fig. 7. — Distribution of Aff a n / Af r ; se obtained fitting the flares with power-law and exponential models. These ratios 
are independent of redshift. 
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Fig. 8. — (left): Distribution of observing times. The gap at log t ~ 3.5 is due to observing constraints (end of the first 
orbit), (right): Time resolution (BT) as a function of the time since the BAT trigger. The solid curve corresponds to 
BT/f = 0.1 and lies above the large majority of the data points. 
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Fig. 9. — The results of the simulation in the (t, At) plane: here we plot the contours for which we have the same 
detection probability. Black points are the real data, based on Gaussian widths and peaks. The dashed lines correspond 
to the three levels Af = 1, 0.1, and 2 x 10 3 . 
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Fig. 10. — Ratio between the peaks of two successive events for both classes: X-ray flares (crosses) and gamma-ray 
pulses (circles), as a function of the quiescent time between the two events. Solid line shows the mean value, -0.157, 
when the two points with the lowest ratio are ignored; dashed lines show the + lcr region. The outliers are GRB 050724 
(log peak i+( . /peak,. ~ 10~ 2 ) and GRB 050502B (log peak i+( . /peak,. ~ 10~ 4 ). 



34 



16 
14 
12 
10 
8 
6 
4 
2 




X-ray flares 
y-ray pulses 
jpulses+flares 



-2 -1 
Log(peak i+1 /peak|) 







Fig. 11. — Distribution of the ratio between the peaks of two successive events: X-ray flares (red), gamma-ray pulses 
(green), both classes (blue). 
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Time since BAT trigger (s) 

Fig. 12. — Gaussian peaks of the flares as a function of time. The solid line is the best fit, while the dashed lines 
correspond to 95% confidence limits. 
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Fig. 14.— AfMi/Af rise vs r 90 . 
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3 




Fig. 15. — Scatter plot of AF IF- At It values calculated on our flare samples on the kinematically allowed regions 
for afterglow variabilities according to lloka et al.l (120051) . Data are drawn from Tables |2]and[3j we used the FWHM 
of the Gaussians as At, the Gaussian peak time for f, while the ratio of the peak flux over the underlying continuum 



flux (A F/F) was calculated using the best fit models. The four limits plotted are based on (a) eq. (4) in lloka et al. 



(2005) for dips (shown on axis), (b) eq. (7) for bumps due to density fluctuations (on axis), (c) At > t for bumps due 



to patchy shells, and (d) Af > f/4 for bumps due to refreshed shocks. According to lloka et alj (12005b . when many 
regions fluctuate simultaneously, limits a and b are replaced by eqs. (Al) and (A2) in Ioka et alT ( 2005b . respectively. 



The off-axis cases (viewing angle 6 V 
shown by dashed lines. 



y /2 > A9, where AO > y is the half-angular size of the variable region) are 
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